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Abstract 

A procedure for reducing the functional integral of QED to an integral over 
bosonic gauge invariant fields is presented. Next, a certain averaging method for 
this integral, giving a tractable effective quantum field theory, is proposed. Finally, 
the current-current propagator and the chiral anomaly are calculated within this 
new formulation. These results are part of our programme of analyzing gauge 
theories with fermions in terms of local gauge invariants. 



1 Introduction 



This paper is part of our programme of analyzing gauge theories in terms of physical 
observables (i.e. gauge invariants). For applications of this programme to non-Abelian 
Higgs models we refer to [ij]. In recent years we have applied it to theories of gauge 
fields interacting with fermionic matter fields, see - In we have proved that 
the classical Dirac-Maxwell system can be formulated in a spin-rotation covariant way in 
terms of gauge invariant quantities. In we have shown that similar constructions work 
on the level of the (formal) functional integral of QED, where fermion fields are treated 
as anticommuting (Berezian) quantities, and in we have applied our procedure to the 
2-dimensional Schwinger model. As a result we obtained a functional integral completely 
reformulated in terms of local gauge invariant quantities, which differs essentially from 
the effective functional integral obtained via the Faddeev- Popov procedure [Q. 
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The present paper is a continuation of 0. It turns out that our general procedure 
leads to a complicated, singular functional integral kernel. In order to make this model 
tractable, we propose a certain averaging procedure leading to an effective quantum field 
theory. This theory is characterized by a certain number of parameters, which have 
to be fixed by comparison with experimental data. As an application we discuss - for 
the massless case - the current-current propagator and the chiral anomaly within this 
formulation. (In principle, the massive case can be also dealt with, using an expansion 
in the mass parameter, but this problem will be not addressed in this paper.) A number 
of interesting phenomena and results comes out: A bosonization rule, similar to that in 
the Schwinger model appears naturally. Moreover, we get a dynamical mass generation 
leading to a massive spin-1 field. Physical quantities like the current -current propagator 
and the chiral anomaly are given as expectation values with respect to an effective non- 
local measure. This measure can be analyzed in terms of a power series expansion in the 
coupling constant, which, however, is completely different from the ordinary perturbation 
expansion. This is due to the fact that the above mentioned mass itself contains the bare 
coupling constant. Therefore, the formulae obtained suggest that automatically some 
resumation of the ordinary perturbation series has taken place. It is also remarkable 
that our formulation leads to a completely new approach to calculate the chiral (Adler- 
Bardeen) anomaly. In lowest order this quantity can be calculated analytically. Adjusting 
part of our free parameters yields the standard Adler-Bardeen anomaly with the correct 
coefficient. 

We stress, that our approach circumvents any gauge fixing and, therefore, also the 
Gribov problem [Q. It leads naturally to bosonization and can be viewed as a general 
construction scheme for effective quantum field theories. Due to the above remarks, it 
seems to be appropriate for the study of non-perturbative aspects. We also mention that 
a similar construction is possible for lattice models within the Hamiltonian framework. 
In this context we have discussed the charge superselection structure of QED [0. 

Our method applies also to Yang-Mills theories, see |T^, where the functional integral 



of one-flavour chromodynamics is reduced to an integral over purely bosonic invariants. 



2 QED in Terms of Local Gauge Invariants 

The functional integral of QED is given by 

Z = y JJdAdV^d^*e*^^'^'^[^''^''^*l, (2.1) 

= -^F^^F^^'' -mr* Pabi^' - Irn[r* Pabhn\D,r} , (2.2) 

where F^^^, = d^Ay — d^A^ and D^'tj)"' = d^i/j'^ + ieA^ip"" are the electromagnetic field 
strength and the covariant derivative, respectively. Here a,b, ... = 1, 2, 1, 2 denote bispinor 
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indices and /U, z^, ... = 0,1,2,3 spacetime indices, (3ab denotes the Hermitean metric in 
bispinor space and (7^^)^^ Dirac matrices. The anticommuting components of the 

bispinor field "0", which can be represented by a pair of Weyl spinors = I 

generate a Grassmann-algebra of pointwise real dimension 8. The representation used for 
these quantities can be found in p| and P]. 

In 0] we have proposed a procedure which reduces the functional integral ( |2.1| ) to 
an integral over gauge invariants. It is based upon the following ideas: First one has to 
analyse the algebra of Grassmann-algebra-valued gauge invariants, which can be built 
from the gauge potential and the anticommuting matter fields Typically, there 
occurs a number of identities between the invariants which, in general, cannot be solved on 
the algebraic level. In particular, one finds a relation which expresses the Lagrangian, or 
the Lagrangian multiplied by some non-vanishing element of the above algebra, in terms 
of invariants. In a next step one has to implement this relation under the functional 
integral and to reduce the original functional integral measure to a measure in terms of 
invariants. For that purpose we make use of the following notion of the ^-distribution on 
superspace (see [§] and [§]) 

^^^2niau-u) ^ y- LJ^5W(^i)[/". (2.3) 

n=0 

Here u is a. c-number variable and U an element of the Grassmann-algebra built from 
matter fields -0 and ip*. iFiom this definition we get immediately 

1 = / du6{u-U). (2.4) 



Thus, by inserting identites of the form ( |2.4| ) under the functional integral, we introduce 
for each Grassmann-algebra-valued gauge invariant a c-number variable, which we call 
c-number mate. These mates are by definition gauge invariant. This way we are able to 
solve the above-mentioned relations between invariants, leading to a theory reformulated 



in terms of gauge invariant fields. For details we refer to B and |10 



For QED we start with the following Grassmann-algebra-valued gauge invariants: 
H 



J" 

•^5 



Vk<P'', (2.5) 
Im {H* {^*kD,4>'' + ^""D,^*^) ] , (2.6) 

r*f^a,m\r (2.7) 



Here is a complex scalar field, whereas is a real-valued covector field. and 
denote the vector and axial- vector current, respectively. 
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We denote the corresponding c- number mates by h, 6^, and j^, and put := 2e\h\^ 
as well as /i = |/i|e*°. It was shown in P] that our procedure yields the following functional 
integral: 

Z = [ ]^{dz;^dj'^dj5^d|/ipdair[i^,j5^J/i|2]}e^^<^'"^l'''^'^''^'5'l^^ (2.9) 

1 f 52 52 ^2 ^2 52 ^2 

J \- 2 4 

16 vr y.dji'dj^dfdj^ ^j'^^U^H^il ^j^^il^f^H 
^2 ^2 1 ^4 1 ^3 1 ^ 

1 5 _1 6^ _ ^ ^ 

^16\h\^ S\h\~2 S\h\^6j^'6J^, ~ 2\h\ 6\h\ Sjt'Sj^ 

C[v,,f,j^,\hla] = -l(di,v,^f-efv, + lj^id,a) (2.11) 

+ she''^^^ j^j^id^j,) -2m\h\cosa. 



where 



and 



3 Effective Bosonized QED 

Observe that the integral kernel K[j^'J^, \h\^] has the form K = V {S{jf)d{j^)d{\h\'^)}, 
where D is a differential operator containing functional derivatives with respect to 
and \h\, multiplied by singular coefficients. A priori, this expression does not make sense. 
In order to regularize it, we replace it by a Gaussian measure with three free parameters, 
aj, ajB and ah, which have to be fixed by physical requirements later. This way we get: 

Z = ^^ /" JJ{dt;^d/di5^d|/i|2da}e'-^^'"^["'"^'''^5M'*l--l, (3.12) 

where now 

-2m \h\ cosa - ^^fj, - ^jt^j^ - jt;Jh\' . (3.13) 

It is interesting to note that this regularization can be achieved by a technical trick 
similar to that used in the Faddeev-Popov procedure: One can average the singular kernel 
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( 2.101 ) with a functional depending on three auxihary fields corresponding to the variables 
j^ii'^ and \h\. It was shown in [|l^ that the requirement to obtain the above Gaussian 
measure after this averaging determines this functional uniquely. 

We see that the Lagrangian ( p.l3| ) does not contain derivatives of the chiral current 
j^, i.e. enters the theory as a non-dynamical field. Thus, we can carry out the simple 
Gaussian integral over the chiral current, which yields 

Z = M y JJ{dv^di'^d|/i|2da}e'^<^'^^['''"^'''l''l'"l , (3.14) 



where the effective Lagrangian has taken the form: 



a,5 



' 128FF ^^MP<x Jaidph)j\d^r) (3.15) 

Comparing this Lagrangian with ( p.l3| ) we observe that has been replaced by the 
gradient of a. More precisely, we have the remarkable relation 

d,a = —jl. (3.16) 

This is the 4-dimensional analog of the bosonization rule in the 2-dimensional Schwinger 
model, see [0. (In |^ we have obtained this rule for the Schwinger model using our 
approach.) 

Observe that the field l/ip enters (|3.15|) in a non-dynamical way, too. Thus, in prin- 
ciple, one should integrate it out. This, however, can not be done explicitly. In a first 
approximation, could be replaced by a constant. We will come to that point later. 

Finally, let us write down the generating functional integral of our effective theory 

dm), (my- 

Z[C^,e^,r7^] =Af f JJ{dt;^di^d|/i|2da}e'^'^'"^[^'''«''''''''"'"^'''l''''"l , (3.17) 

with 

+^ {d^a){d^a) - 2m\h\ cosa + C{d^a) + ^^j^ + r]^v^ , (3.18) 
where and rj^ denote the source currents for (d^a), and v^^, respectively. 
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4 The Current— Current Propagator 

In this Section we want to calculate the current-current propagator < 0\Tj^{yi) j'^{y2)\0 > 
using the effective theory obtained in the last Section. We restrict ourselves to the massless 
case, i.e. we put m = in the generating functional integral ( p.l7| ). 
The current-current propagator is given by 



<0\Tf{y^)f{y2)\0> 

1 1 S 



6 



z[o,o,o] 12 6Uyi) ^Uy2\ 

115 5 



z[o,o,o] 12 5Uyi) ^Uy2) 

To handle the non-linear (self-interaction) term 

occuring in ( |3.18| ) we introduce the vector field 

and decompose it into its longitudinal and transversal parts: 

k^[J,^h\] = ^^k[3,^h\] + k^^[3,,\h\] 

where d^k^^-^[j^, \ h\] = 0. This yields 



(4.19) 



(4.20) 



(4.21) 



C[C>^,e,v'';v^,f,\h\,a] 



-I (A, 



ej% + {d,a) {d^k[j,, M) + g {d^k[j„ \h\]){d,k[j„ \h\]) 
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k^^[j„ \h\] kilj,, \h\] - ^fj, -^Jh\^ + ^ {d^a){d, 

The term ^ {d^a) k^^[j^, \h\] vanishes by partial integration. 
Transforming 

"' = " + 4^^b'/^'l^l]' 



(4.22) 



(4.23) 



we obtain 



c[o,e,v^;v^,f,\h\,a'] 

= -i id^,vu]? - efv, + 31 \h\] kilj,, \h\] 



(4.24) 



Observe that the transformation ( 4.23 ) leaves the integral measure da invariant. Thus, 
a' can be integrated out trivially and we get 



<oim?/i)/(?/2)io> 
1 



Z[0,0,0] 



n{dt;^dj'^d|/ip} 



1 5 5 



(4.25) 



where 



= -\ {d[,vui)' - efv, + g k^^[j„ \h\] k^ij,, \h\] 

To" "''^■^'^ 2 Oft I 



(4.26) 



The further analysis of the term k^^^[jf^, \h\] ^^[j'/i, \h\] needs some care, due to its 
non-linear (and non-local) character. Expanding the exponential of this term in a series, 
we obtain 



<oimt/i)i''(!,2)|o> 



X 



n=0 

1 ^ jfd'^xCo\0,(.t',v'';Vu,j'',M 



Z[0 



^/nKd/d|/.r 



X 



ra=0 

1 5 

X ^ 



,i/d4x£o[0,C'',r;'';i,^,if,|h|] 



where 



^[0,0,0] = y n{d^Md^''d|^l'E^ (S)" (^"^[iirJ/^P^irlTr,!/^!])" 



X e 



x n=0 
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and 



Co[o,e,v'';v„f,\h\] 



efv^ 



2 a, •^'^■^M 



2a 



In the last step we used the fact that the integral kernel consists of a polynomial function 
in j'^. Thus, j'^ can be replaced by the corresponding functional derivative with respect 
to the source current Now we are left with a Gaussian integral with respect to j^, 
which yields 

<0|T/(?/i)/(l/2)|0> 

-IP °° 

' ' J f n=0 



Z[0, 



X 



6 



6 



J J d4a;£o[0,C'^,r)'';Dp,|/i|] 



(4.27) 



with 



/oo 

n=0 



X e 



/ d-*x£o[0,e'^,r;'^;i)M,|ft|] 



and 



/:o[o,e^^^;^^,|/i|] 



(4.28) 

We remark, that the covector field has acquired a mass ml = ajc^. Thus, in 
our effective field theory, the original gauge potential has been replaced by a mas- 
sive spin-1 field v^. Now, observe that due to (|4.28| ) the functional derivative with re- 
spect to ^'^ produces the term —ajcv^. Therefore, the non-linear (and non-local) term 

{k^'^ii^j 1^1] ^/i'[r^' 1^1]) is effectively of the order e'^, and we can treat it as a per- 
turbation. Performing the functional derivatives in (|4.27|) with respect to yields a 
complicated non-local measure, which in full generality cannot be handled analytically. 
Limiting ourselves to lowest order we get: 

<0\Tny,)f{y2)\0>o 
1 

- m 



l[{dv,d\h\'} 

X 



(4.29) 
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Here we have denoted 

:C(i[0;v^„\h\] 



^[0] = j Y[{dv^d\h\^]e'^^"-' 



and 

CoV)^- v„ \h\]^-\ {d[,v,])' + 4 ^;% - ^Jh\' + ^% ■ (4.30) 
Next, observe that the integration over \h\ decouples, giving 

<o|r/(i/i)/(?/2)|o>o 

= < aj rj"" 6\yi - ^2) + mlaj v''{yi)v''{y2) > , 
where now 

Z[0]= / J]{d^;4e'/^'^^°[°'''''l 

and 

This way we obtain: 

Proposition 1. In lowest order, the current-current propagator < OlTj'^^x) j'^{y)\0 > of 
massless QED is given by the vacuum expectation value 

< 0\Tnyi)r{y2)\0 >o=< cXjV'"'S\yi-y2)+mlajv'^{y,y{y2) > (4.31) 

with respect to the functional measure Yix {d'f//} e^-f '^^^^^^'^^ and the Lagrangian 

C[v,\^-\{d^,v^^f + 4v^v,. (4.32) 
Thus, in lowest order, we are left with a simple Gaussian integration. We have 

with D""^ :— d'^dfj, — S'^^d'^d^ — ml 5^. We introduce the new field 

v'^'ix) - v^(x) + J d'yr)''(y)(D-'r^(x-y) , 

where the propagator {D~^)^^^{x — y) of the massive free field v^j, is defined by 

D\{D-'r^{x-y)^-5\x-v)5>i . 



Thus 



iD-'r.ix-y) = ',;^l^6\x-y). (4.33) 

Performing the above transformation, the resuhing integration over v'^ decouples and we 
obtain 



<0|T/(2/i)/(l/2)|0> 







%(?/i) 5??v(2/2 

lr;''=0' ^ ' ' 

Finally, performing the remaining functional derivatives, we get 

<0|r/(2yi)/(2/2)|0>o = a,^^^ b\y,-y,)-a,'-^^^^^^^\y,-y,) 

= c.,^-^^^6\y,-y,). (4.35) 

Fourier transforming to momentum space leads to the following expression in lowest order: 

nr(p) = < OlTfiyi) ny2)\0 >o= {p'p' - v'"?') Tip') (4.36) 

with T{p^) := ^2lp2 ■ This result has the expected Lorentz structure. Moreover, we obtain 
the identity IIq'^{p) = 0, which is nothing but the vector Ward identity. Thus, in lowest 
order of the above defined perturbation series, our result obeys the vector Ward identity. 

However, a direct comparison with the well-known perturbation series of QED is not 
possible. This is due to the fact that the mass ml = a^e^ of the spin-1 field occurring in 
( 4.35| ) contains the bare coupling constant e. Thus, expanding this result around = 0, 
which corresponds to a power expansion in e^, we obtain non-vanishing contributions 
to all orders in e^. Therefore, formula ([4.35|) can be interpreted as a resumation of 
particular quantum corrections, which results in an effective ("dynamical") mass for the 
field Vfj,. Unfortunately, higher order contributions (in the sense of our expansion) cannot 
be calculated analytically. 



5 The Chiral Anomaly 



In this section we show that within our approach the correct Adler-Bardeen anomaly |n 
is obtained in the lowest order approximation discussed in the previous section. Again we 
restrict ourselves to the massless case. 

Due to the bosonization rule ( |3.16| ) the chiral anomaly < d^j^ > is given by the 
vacuum expectation value 

<d,j^> ^ <^{d,d^a{y))> 

^ T (^T7^) Z[C,e,v'] , (5.37) 



Z[0,0,0] 2i V '^SCf^iy) 
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where ^f^, 77^] is given by ( ^.17| ) and ClC^, 77^; t;^, j^, \h\,a] by ( |3.18| ), or in terms of 
k'^ljfj,, \h\], by ( [4.22|) . Choosing a longitudinal source current C'^ := {d^O and transforming 

«" = « + ^C+4^MJm,I^I] (5.38) 
yields the Lagrangian in the following form: 

= - efv, \h\] + S^'^b.' Mkilj,, M 

Now a" can be integrated out trivially. Performing the functional derivative with respect 
to C'^ we obtain 

a,-5 



< ^ {d,d^a{y)) > 



Z 



X e 



i / d4x£[C^,^'',0;i)pJ'',|h|] 



X e 



/d4x£[0,5'^,0;t;^,j'^,|h|] 



^M=0 



where 

C[Q,e,0;v,,f,\h\] 

= -\ {dy,v,\f - efv, + 31 k'^^ij,, \h\]ki[j„ \h\]-^ fj, -^Jh\^ + ej.. 

Due to the explicit dependence of k'^[j^{y),\h{y)\] on \h\, a further exact analytical 
treatment of this formula is impossible. But, as outlined in Section ^, \h\ enters our 
effective theory in a non-dynamical way and, therefore, in principle it can be "averaged" 
out. In the simplest approximation we replace \h\ by its mean value Xo- This way we are 
led to 



< ^ {d,d^a{y)) > 
1 



Z[0,0,0] 



l[{dv,df} {-ld,k^[j,{y),Xo]) 



X e 



Z[0,0,0] 



l[{dv,df}(-^e'^^^{d,js{y)mjM) 



X e 



:f=0 



(5.39) 
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with 



Now the non-trivial coupUng term '^k^^'^[i^^xo\ ^/[[j/iiXo] will be treated similarly 
as in the previous Section. In lowest order we get 

< ^ {d,d^a{y)) >o 



Z[0 



X e 



-14 



i/d4x£o[0,«'',0;^;^,j''] 



where 



Z[0, 0, 0] = /" JJ jdv/, dj'^} e'/'i'^^otO'O'O'"'' 



and 



Performing the Gaussian integration over and the functional derivatives with respect 
to we finally obtain 



8xg 



2^2 



e 



<(a[^t;,](2/))(a[„t;^] (!/))> (5.40) 



with Z[0] = /rix {d^^} e^/'i'^^l'^^l and 

>C[^M] = -i(AM^.])' + ^^%- (5.41) 

Thus, bearing in mind that F^j^^, = Sj^fi^] we can formulate the following result: 

Proposition 2. In lowest order, the chiral anomaly of massless QED in (3 + 1) dimen- 
sions is given by the vacuum expectation value 

< ^ (d^d'^aiy)) >o= ^ e^''^^ < F,, F^p > (5.42) 
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with respect to the functional measure 

nxid^^Mje'-^'**''^^''^' , where 



(5.43) 



If we start with an external electromagnetic field, becomes external, too. In that 
case we get 



jivafi rp p 



(5.44) 



where F^^, = d[^Vy-^ denotes the external electromagnetic field strength. 

Observe that we get the correct coefficient, see [|13[, |T^, for the chiral anomaly if we 
require the following relation between our parameters: 



Xo 



1 

2 7r2 • 



This yields 



<d^jl{y))>.^<^{d,d^^a{y))> 



16 TT" 



C -L li.lj ^ ( 



(5.45) 



(5.46) 



Thus, the anomaly can be used to fix one of the parameters of our effective theory as 
discussed in the Introduction. 
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